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摘 要 : 为 提高 Hargreaves-Samani( H-S) 模型 计 


下 游 平 原 共 128 个 气象 站 点 1961—2010 4 hy 2 


参考 蒸 散 发 的 精度 ,利用 西北 黄河 流域 与 长 江 中 
日 气象 资料 对 H-S 模 型 进行 回归 修正 ,以 Penman- 


Monteith(P-M) 模 型 为 标准 ,评价 了 H-S 改 进 模型 H-Scom 模 型 的 计算 精度 ,并且 以 第 六 次 国际 耦合 
模式 比较 计划 (CMIP6) 气 候 模式 来 对 HScom 模 型 进行 了 未 来 适应 性 评价 。 结 果 表 明 : 修 正 后 ,在 
验证 期 内 ,长 江 中 下 游 平 原 4 个 分 区 的 平均 绝对 误差 (MAE) 和 和 均 方 根 误差 (RMSE ) 的 平均 值 分 别 下 
降 了 6.21 mm: AF 6.38 mm 月 "1; 西北 黄河 流域 4 个 分 区 的 MAE 和 和 RMSE 的 平均 值 分 别 下 降 了 
9.26 mm 月 "和 9.23 mm ' 月 ,2 个 研究 区 域 修正 后 的 决定 系数 ( 尼 ) 比 修正 前 最 少 提高 1%。 在 
CMIP6 气 候 模 式 的 未 来 气候 情景 下 屁 均 达到 0.98 以 上 ,具有 和 良好 的 适应 性 。 该 研究 修正 的 模型 方 
法 可 为 仅 有 气温 数据 的 地 区 提供 较 高 精度 的 参考 燕 散 发 估算 方法 ,为 高 频 灌 溉 提供 较为 准确 的 数 


据 基础 。 


关 键 词 : 参考 作物 蒸 散 量 ; Hargreaves-Samani 模型 ; 
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参考 蒸 散 发 (ETo 是 气候 \ 水文 ,水 资源 规划 和 
管理 等 诸多 应 用 中 的 一 个 重要 参数 " ,其 准确 估算 
对 于 灌溉 系统 设计 、 水 资源 管理 .农业 和 水 文 气象 
研究 以 及 水 资源 收文 的 确定 至 关 重 要 ”人 。 常 用 的 
ET 计算 方法 有 Penman-Monteith (P-M) 模 型 5 、 
Makkink 法 、Hargreaves-Samani(H-S) 模 型 ”。 目 
前 ,联合国 粮农 组 织 1998 年 提出 的 P-M 模 型 是 较为 
精确 的 ET, 计 算 模 型 人 中。P-M 模 型 被 认为 是 标准 方 
法 ,与 蒸 渗 仪 在 各 种 气候 条 件 下 的 实测 值 和 许多 其 
他 模式 的 实测 值 相 比 ,P-M 模 型 可 以 应 用 于 各 种 环 
境 和 气候 情景 ”。 然 而 ,该 模型 的 计算 过 程 相 对 复 
林 , 需 要 用 到 大 量 气 象 数 据 ”。 包 括 人 研究 尺度 内 最 
高 气温 .日照 时 数 .最 低 气 温 、 风 速 和 相对 湿度 ,这 
限制 了 它 在 许多 地 方 的 应 用 ,因此 需要 基于 气温 等 
少数 气象 要 素 的 ETo 计 算 方 法 。H-S 模 型 是 目前 应 
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用 最 广泛 的 基于 温度 的 方法 之 一 ,基于 H-S 模 型 
对 ET 的 准确 估计 ,可 以 实现 ET 的 短期 预测 H-S 
模型 虽然 简单 可 靠 , 但 H-S 模 型 的 准确 性 往往 受到 
当地 气象 数据 的 影响 ,不 仅 是 最 高 和 最 低温 度 ,还 
有 风速 和 湿度 ”。 研 究 发 现 ,在 不 同 的 气候 区 域 ， 
H-S 模 型 的 校正 系数 每 月 都 有 差异 。 而 实际 应 用 
中 ,部 分 地 区 观测 的 气象 数据 又 很 难 满 足 P-M 公 式 
的 数据 要 求 "“"。 因 此 ,结果 精确 、 适 用 范围 广 的 ET。 
简化 计算 方法 得 以 发 展 。 对 H-S 模 型 进行 修正 即 可 
以 提高 计算 精度 ,对 未 来 的 ET 估算 也 有 积极 意义 。 
迄今 为 止 ,在 全 球 不 同 气 候 条 件 下 ,H-S 模 型 的 
校准 和 验证 已 经 开展 了 多 项 研究 。Gavilan 等 " 根 
据 瓜 达尔 基 维 尔 山谷 灌溉 区 的 数据 证 明 ,如 果 使 用 
修正 后 的 方程 而 不 是 原来 的 方程 ,用 水 量 可 以 减少 
11%。 因 此 ,为 了 更 好 地 适应 局 部 区 域 ,有 必要 对 
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H-S 模 型 进行 校准 。 冯 克 鹏 等 “通过 最 小 二 乘 回归 
对 H-S 模 型 进行 修正 ,确定 出 宁夏 地 区 每 个 气象 站 
点 的 经 验 系数 ,从 而 得 出 修正 后 的 计算 公式 。 夏 兴 
生 等 ”采用 线性 回归 分 析 方 法 对 H-S 模 型 的 系数 进 
行 逐 月 逐 区 的 修正 ,并 以 最 优选 系数 确定 了 不 同 分 
区 的 模型 校正 系数 。 这 些 工 作 只 显示 了 参数 和 性 
能 指标 值 分 布 的 一 般 规 律 ,缺乏 日 稼 评价 和 大 量 的 
站 点 。 为 了 在 高 频 灌 溉 中 应 用 更 精确 的 ETu,H-S 模 
型 需要 使 用 来 自 全 国 众多 站 点 较 长 的 日 气象 数据 
序列 进行 校准 和 验证 。 

本 人 研究 通过 长 江 中 下 游 平原 与 西北 黄河 流域 
共 128 个 气象 站 点 1961 一 2020 年 的 逐日 数据 对 H-S 
模型 进行 修正 ,提高 该 模型 对 气候 因子 随时 间 变 化 
的 响应 能 力 。 并 深入 分 析 比 较 了 该 模型 在 西北 干 
旱 半 干旱 区 (西北 黄河 流域 ) 和 湿润 区 (长 江 中 下 游 
平原 ) 的 计算 精度 与 适应 性 ;并 对 校准 后 的 H-S 模 型 
进行 未 来 的 适应 性 分 析 。 


1 研究 区 概况 


长 江 中 下 游 平 原 主要 由 江汉 平原 洞庭 湖 A 
阳 湖 、 皖 苏 沿江 、 里 下 河 及 长 江 三 角 洲 6 块 平原 组 
成 ,地 跨 湖北 、 湖 南 、 江 苏 、 江 西 浙江、 安徽、 上 海 7 
省 市 (图 1) ,海拔 在 5~100m 之 间 。 其 大 部 分 位 于 北 
半球 的 亚热带 ,多 年 降水 量 为 1000~1500 mm, 多 年 
平均 气温 14~18 % ,多 年 平均 ET, 为 1040.91~ 
1091.65 mm。 该 地 区 较 低 的 ET 和 较 多 的 降雨 量 使 
其 气候 较 同 纬度 内 陆地 区 更 为 湿润 , 亦 为 中 国 水 资 
源 最 丰富 的 地 区 。 该 地 区 降雨 量 和 ET 季节 分 布 较 


(a) 长 江 中 下 游 平原 


0 80km 


均匀 ,日 与 国内 大 部 分 湿润 区 的 ET, 和 降雨 空间 分 
布 相似 。 

西北 黄河 流域 的 地 势 西 高 东 低 ,西部 河源 地 区 
平均 海拔 4000 m 以 上 ,地 形 地 貌 以 高 原 为 主 ,冰川 
积 雪 遍布 。 其 大 部 分 地 区 年 降水 量 在 200~650 mm 
之 间 , 年 平均 气温 为 12~14 CC。 夏季 气温 高 ,降水 稀 
少 , 且 自 东南 向 西北 逐渐 递减 。 冬 季 寒 冷 干燥 ,水 
资源 相对 匮乏 ,生态 环境 脆弱 。 该 地 区 多 年 平均 
ET) X 845.93~1225.07 mm , 较 其 他 半 干 旱 区 ,该 研究 
区 域内 流域 蒸发 能 力 较 强 ,甘肃 .宁夏 地 区 属国 内 
年 蒸发 量 最 大 的 地 区 ,是 中 亚 干 旱 区 的 重要 组 成 部 
分 , 占 国内 干旱 半 干 旱 区 的 60%。 


2 数据 与 方法 


2.1 数据 来 源 

本 人 研究 利用 中 国 气象 资料 服务 中 心 (http://data. 
cma.cn) 提 供 的 西北 黄河 流域 与 长 江 中 下 游 平原 共 
128 个 气象 站 点 1961 一 2020 年 的 最 高 气温 、2 m 高 
处 风速 .降雨 量 日照 时 数 、 相 对 湿度 和 最 低 气 温 等 
气象 资料 (站 点 位 置 如 图 1 所 示 )。 对 所 有 站 点 ,分 
别 使 用 1961 一 2010 年 .1981 一 2010 年 的 数据 进行 回 
归 修 正 ,2011 一 2020 年 的 数据 进行 验证 。 

降雨 是 水 输入 生态 系统 的 主要 方式 之 一 ,水 分 
条 件 影响 能 量 成 分 之 间 的 净 辐 射 分 布 比例 ,例如 显 
热 通 量 和 潜 热 通 量 ,从 而 改变 空气 水 热 条 件 和 水 气 
FR BE GR , SPL XT ET WSR AD PEO, AN CAR BE 
1961 一 2020 年 128 个 气象 站 点 的 年 均 降 雨量 数据 ， 
采用 反 距离 空间 插值 对 人 研究 区 域 进 行 降 雨量 插值 ， 


(b) 西北 黄河 流域 


>z 


注 : 该 图 基于 国家 测绘 地 理 信息 局 标准 地 图 服务 网 站 下 载 的 审 图 号 为 GS(2019)1822 号 的 标准 地 图 制作 , 底 图 边界 无 修改 。 下 同 。 


图 1 气象 站 点 位 置 示意 图 


Fig. 1 Schematic diagram of the location of meteorological stations 
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分 别 对 2 个 研究 区 范围 内 的 降雨 量 进行 均等 分 区 (图 
2) ,以 讨论 不 同 水 分 条 件 区 域 ETo 计 算 精 度 。 
2.2 H-S 模型 的 校准 
2.2.1 参考 蒜 散 发 计算 方法 

(1) P-M 模型 

P-M 模 型 被 FAO 推 荐 为 计算 ET 和 评价 其 他 
方法 的 标准 方法 。 因 此 ,将 观测 气象 数据 输入 模 
型 ,利用 输出 的 ETo 结 果 量 化 H-S 模 型 计算 的 ETo 精 
度 。P-M 模 型 计算 ETo( ET,_py ) 采 用 FAO 官网 推荐 
的 ETo Calculator 软件 ,其 计算 公式 为 : 
0.408A(R, - 6) + 20 _yu,(e,-e,) 

A+y(1+0.34U,) 

RP: Ry sed ae M md);C 为 土壤 热 通 
(MJ md) ;7 为 平均 气温 (%C);y 为 湿度 计 和 常数 
(kPa ;号 为 2m 高 处 风速 (ms ie. 为 饱和 水 
UK (kPa); e, 为 实际 水 汽 压 (kPa); A 为 饱和 水 汽 
压 与 温度 曲线 斜率 (kPa:C7)。 

(2) H-S 模 型 

H-S 模 型 是 Hargreaves 等 "在 198$ 年 提出 的 经 
验 公 式 , 与 P-M 模 型 相 比 ,H-S 模 型 计算 更 简单 ,对 
气象 数据 的 要 求 更 少 ,因此 该 方法 可 以 在 世界 范围 
内 广泛 应 用 。H-S 模 型 计算 ETo ETo- ys ) 公 式 如 下 : 

ET 135 = 0.00238, (Ty, + 17.8)(T a- Tan) (2) 


ave 


ET og = (1) 


R,= Ax 6 d,(o, sin $ sind + cos b cos Ôsin w,) (3) 


AP: R, 为 天 顶 辐射 (MJ md’); Thos Tan > Tove 
分 别 为 最 高 温度 .最 低温 度 .平均 温度 (%); Go 为 太 


(a) 长 江 中 下 游 平原 


>z 


T 9 soim 
降雨 量 /mm 
mm A 区 (1640~1919 mm) 
EE BK(1361~1640 mm) 


E C 区 (1082~1361 mm) 
OO DK(804~1082 mm) 


KA Fe BL , AS SCHL 0.082 MJ.m .min'; d 为 日 地 相对 
PERS (km); w. 为 日落 时 角 (°);6 为 太阳 倾角 (°), 与 
每 天 在 一 年 中 的 序数 J 有 关 ; 4 为 纬度 所 转换 的 弧 
度 (rad)。 其 中 ,6、d, Q o, 等 参数 按照 下 式 进行 
HA: 


z in 27 1- 
5=0.409 sin( 2E J 1.36) (4) 
5 2m 
d, = 1+0.033 cos 27.) (5) 
= TR0(9) (6) 
w = arccos(—tan ¢ tan6) (7) 


式 中 :7 为 一 年 中 的 时 间 序 数 ;e 为 当地 纬度 (")。 
2.2.2 修正 方法 H-S 模 型 相 比 P-M 模型 对 气象 数 
据 的 要 求 更 少 , 且 计算 过 程 简单 ,可 以 更 广泛 的 运 
用 。 但 相对 于 P-M 模 型 ,H-S 模 型 对 ET, 的 计算 精度 
较 低 ,因此 需要 对 H-S 模 型 进行 修正 ,以 满足 精度 要 
求 , 可 以 较 好 的 适用 于 气象 数据 缺乏 的 地 区 。 

基于 非 线性 最 小 二 乘 方法 ,提出 了 对 原始 H-S 
模型 修正 后 的 阳 - Sconn 模型 。 研 究 过 程 主 要 包括 
ET o-ou 和 了 Im 的 计算 (ETw 和 了 Ta 分 别 为 
P-M 模 型 和 H-S 模 型 所 计算 的 ETo) ,H-S 模 型 校正 系 
数 的 计算 以 及 误差 分 析 , 以 及 对 2021 一 2030 年 的 
ET 预测 。 本 研究 选择 1961 一 2010 年 和 1981 一 2010 
年 2 个 校准 周期 作为 校准 期 ,2011 一 2020 年 为 验 
证 期 。 

以 各 气象 站 点 多 年 的 逐 月 ET,» 为 因 变 量 ， 
ET, 4; YASH. 建立 修正 方程 ET sy =a+b 


(b) 西北 黄河 流域 


图 例 0 125km 
降雨 量 /mm 

EE EX(700~927 mm) E G 区 (246~473 mm) 
MS FX(473~700 mm) [I HiX(19~246 mm) 


图 2 降雨 区 域 划分 示意 图 


Fig. 2 Schematic diagram ofthe division of rainfall areas 
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ET) ys ,以 此 计算 得 出 矫正 系数 cc ,再 通过 方程 
ET, ys =a+b ET us 结合 2011 一 2020 年 的 气象 资 
料 对 H-Scom 模 型 进行 验证 4 ETo_as ， 为 了 Scom 模 型 
所 计算 的 ET,)。 并 对 比 2021 一 2030 年 H-Scorx 模 型 
与 P-M 模型 计算 的 ETu, 对 H-Scom 模 型 进行 未 来 适 
用 性 评价 ,算法 改进 流程 见 图 3。 


ed 


‘| 196120104 BH 
‘1 1981 一 2010 年 数据 | 


注 :ETo 为 参考 蒸 散发 ;P-M 模 型 为 Penman-Monteith 模型 ;H-S 模 型 
为 Hargreaves-Samani 模 型 ; ET，, 和 ET 分别 为 P-M 模 型 和 H-S 
模型 所 计算 的 ETo; ET yas, A H - Scone! SETHE AY ET); GCM 为 
大 气 环流 模式 ;ACCESS-CM2 为 第 六 次 国际 耦合 模式 比较 计划 
(CCMIP6) 气 候 模式 ;MAE 为 平均 绝对 误差 ;RMSE 为 均 方 根 误差 ; 尽 
为 决定 系数 。 下 同 。 
图 3 算法 改进 流程 图 


Fig.3 Algorithm improvement flowchart 


23 精度 评估 标准 
本 文 以 P-M 模 型 为 参考 标准 ,采用 平均 绝对 误 
差 (MAE) 、 均 方 根 误 差 (RMSE) RE HBR?) 3 S 
§ 标 比较 H-Scomx 模 型 与 P-M 模 型 计算 ET 之 间 的 吻 
合 程度 。 其 中 , 尼 能 够 反映 H-Scom 模 型 与 P-M 模型 
二 者 估计 值 的 相关 程度 ,表明 两 者 的 变化 趋势 和 分 
布 规律 是 否 一 致 。MAE 和 RMSE 能 够 反映 二 者 的 
偏离 程度 。 计 算 公 式 如 下 : 
ET 


2 ET) -us ~ 4" 0- PM(i) 
MAE = = (8) 


n 


Pa ~ Pa) 
RMSE= 


2 


= (9) 


n 


2 2 (iay ~ ET, sci (ET sui 7 ET, om) 
R’ =— 2 


: 2 
py (seen T ET o) (Eia noT BT op vo) 


i=l i=1 


(10) 

AF sn AD REAR Es ET ugg Eag 25 H-S 
模型 计算 的 ET 计算 值 及 平均 值 (mm + H’); 
ET。，wo、ET。awwo 分 别 为 P-M 模 型 计算 的 ETo 计 算 值 
及 平均 值 (mm 月 ");i 为 ETo 的 样本 序号 。 
2.4 H-Sor 模型 未 来 适应 性 评估 

SSP2-4.5 是 第 五 次 国际 耦合 模式 比较 计划 
(CCMIP5 ) 中 RCP4.5 情景 的 更 新 版 本 , 它 结 合 了 一 个 
中 等 社会 脆弱 性 和 中 等 辐射 强迫 的 情景 ,该 情景 还 
稼 被 用 于 区 域 降 尺 度 和 年 代 际 气候 预测 , 故 选 用 此 
场景 ,与 本 研究 更 为 契合 。 以 CMIP6 气 候 模拟 实验 
下 的 ACCESS-CM2 模式 ,以 及 未 来 试验 情景 SSP2- 
4.5 对 HH- Seon 模型 进行 未 来 适应 性 评价 。 对 比 
2021 一 2030 年 H-Scom 模 型 计算 的 ET 与 P-M 模 型 计 
算 的 ET。, 进 行 精度 评估 。 


3 结果 与 分 析 


3.1 ET, 空间 变化 分 布 规律 

从 图 4 可 以 看 出 ,长 江 中 下 游 平原 的 ET 空间 分 
布 呈现 由 南 向 北 递 增 的 趋势 ,为 1041~1092 mm。 其 
中 , 老 河口 和 房 县 2 地 是 ET 的 高 值 中 心 ,为 1081-~ 
1089 mm, 主要 由 于 该 地 地 形 等 地 理 环 境 的 影响 ,对 
于 风速 和 日 照 时 数 的 敏感 性 较 高 ;永州 .常宁 、 道 县 
等 地 低 值 中 心 为 1041~1065 mm, 主要 归 因 于 这 3 个 
地 区 相对 湿度 较 大 日照 时 间 较 短 。 西 北 黄河 流域 
ET 为 846~1225 mm, ,空间 分 布 从 西 癌 东 呈 现 增 加 的 
趋势 。 其 中 惠 农 、 陶 乐 .榆林 、 横 山 和 盐池 等 地 的 多 
年 平均 ET 最 大 ,为 1030.62~1225.07 mm, 这 些 地 区 
日 照 时 间 较 长 ,是 西北 日 照 时 数 的 高 值 中 心 , 且 相 
对 湿度 较 小 ,因此 ETo 较 大 。 
3.2 HH- Soonn 模型 精度 分 析 

原始 H-S 模 型 .H-Scom 模 型 和 P-M 模 型 计算 ET。 
的 统计 指标 (包括 MAE RMSE 、R) 误 差分 析 如 图 5 
所 示 。30 a 校准 期 和 50 a 校准 期 分 别 表 现 出 不 同 的 
精度 误差 ,经 1981 一 2010 年 校准 过 后 的 H-Scons 模 型 
精度 明显 低 于 1961 一 2010 年 校准 模型 所 计算 的 精 
度 , 说 明 更 长 时 间 序 列 数据 的 校准 周期 可 以 减少 气 
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图 4 采用 P-M 模 型 计算 的 ET 空间 分 布 
Fig. 4 Spatial distribution of the ETs calculated by P-M model 


候 变 化 引起 的 误差 和 验证 应 用 存在 的 误差 。 长 时 
间 序 列 数 据 对 于 模型 的 修正 更 加 有 利 ,其 计算 结果 
更 加 具有 稳定 性 。 因 此 ,H-Scom 模 型 采用 1961 一 
2010 年 为 校准 期 。 

经 50 a 校准 期 修正 后 的 H-Scont 模 型 其 MAE、 
RMSE. 尺 在 8 个 分 区 都 有 较 好 的 改善 ,其 相对 应 的 
评价 精度 都 有 明显 的 提高 。 西 北 黄河 流域 的 4 个 分 
区 (RMSE=4.18 mm 月 "') 整 体 比 长 江 中 下 游 的 4 个 
分 区 (RMSE=4.37 mm: 月"') 修 正 效 果 要 好 。 且 修正 
后 的 MAE 值 下 降幅 度 前 者 (9.26 mm 月 ) 比 后 者 
(6.21 mm: A!) EK , WPH H-Scom 模 型 在 西北 黄河 
流域 适应 性 较 强 。 

Be 1 为 原始 H-S 模 型 和 修正 后 H-Scom 模 型 的 统 


(a) MAE 和 RMSE 前 后 变化 


me 修正 前 MAE 上 国 修 正 前 RMSE 
国 修 正 后 MAE (50 a 校准 期 ) 国 修 正 后 RMSE (50 a 校 准 期 ) 
mm 修正 后 MAE (30 a 校准 期 ) m 修正 后 RMSE (30 a 校准 期 ) 


小 
n 


MAE/mm: H, RMSE/mm: H 
= Ww 
T oS 


| | | 


G 


| 


分 区 


出 


| T 


ahh I | | 


0.8 
0.7 
o6 | ™ 修正 前 RR 
i mm 修正 后 R? (50 a 校准 期 ) 
mm 修正 后 R? (30 a 校准 期 ) 


计 指 标 比 较 。MAE 和 RMSE 的 最 小 值 分 别 最 少 降 
低 了 0.60 mm .月 :和 0.56 mm: H~, eK MAE 和 
RMSE 分别 至 少 减少 了 9.00 mm' 月 和 9.28 mm: H“, 
MAE 和 RMSE 平 均值 分 别 下 降 了 7.73 mm :月 和 
7.81 mm* 月 ,修正 后 RR 值 大 部 月 份 均 大 于 0.9,RR 最 
少 提 高 1%。 这 些 结果 表明 ,用 修正 后 的 参数 计算 
出 的 ETo 值 更 为 准确 ,说 明 修正 后 H-Sconr 模 型 可 用 
于 该 研究 区 域 ET 的 估算 。 对 于 西北 黄河 流域 ,经 
修正 后 ,大 部 分 地 区 的 MAE 在 1.11~8.94 mm: HZ 
间 。 对 于 长 江 中 下 游 平 原 ,H-Scom 模 型 的 精度 比 西 
北 黄河 流域 略 低 , 修 正 后 的 效果 略 不 足 ,H-Scom 模 型 
在 西北 黄河 流域 的 精度 较 高 。 总 体 而 言 ,H-Scorr 模 
型 在 研究 区 域内 的 表现 均 明 显 强 于 原始 H-S 模 型 。 


(b) R 前 后 变化 


, | | 


注 :MAE 为 平均 绝对 误差 ;RMSE 为 均 方 根 误差 ; 尺 为 决定 系数 。 
图 5 H-S 模 型 改进 前 后 与 P-M 模 型 的 MAE .RMSE Rz 变 化 
Fig. 3 Changes of MAE, RMSE and R? between monthly ET) and the estimated values of the P-M model 


before and after H-S model improvement 
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表 1 原始 H-S 模 型 和 修正 后 H-Scom 模 型 的 统计 指标 对 比 


Tab. 1 Comparison of statistical indicators between the original H-S model and the revised H-Scorr model 


分 区 ”统计 指标 ”修正 前 后 最 小 值 ” 最 大 值 ”平均 值 


AK MAE 原始 4.08 24.22 12.47 
修正 后 1.76 7.99 3.42 

RMSE 原始 4.31 24.56 13.24 

修正 后 2.23 9.27 4.19 

R 原始 0.71 0.93 0.82 

修正 后 0.90 0.97 0.93 

B 区 MAE 原始 3.18 23.25 11.57 
修正 后 1.53 8.07 3.10 

RMSE 原始 3.78 23.45 12.51 

修正 后 1.96 8.91 3.87 

R 原始 0.71 0.94 0.86 

修正 后 0.83 0.96 0.92 

C 区 MAE 原始 2.50 20.60 8.15 
修正 后 1.69 8.48 4.26 

RMSE 原始 3.22 21.16 9.15 

修正 后 2.08 0.92 5.11 

R 原始 0.66 0.92 0.85 

修正 后 0.82 0.94 0.92 

D 区 MAE 原始 2.74 17.78 7.58 
修正 后 2.14 8.78 4.16 

RMSE 原始 3.56 22.35 8.90 

修正 后 3.00 9.34 5.09 

R 原始 0.65 0.97 0.83 

修正 后 0.81 0.98 0.92 


注 :MAE 为 平均 绝对 误差 (mm 月 ");RMSE 为 均 方 根 误差 (mm: 月 ") 


3.3 H-Scom 模 型 的 应 用 

因 H-Scom 模 型 在 西北 黄河 流域 的 适应 性 较 强 ， 
故 选取 在 西北 黄河 流域 的 下 区 来 进行 了 -Sun 模型 
对 于 2021 一 2030 年 的 ET 模拟 ,以 此 来 验证 H-Scorr 
模型 在 未 来 的 适用 性 。 基 于 CMIP6 模式 数据 评估 
T H-Scom 模 型 在 未 来 的 模拟 能 力 ,并 给 出 了 未 来 
ET。 的 变化 。 本 文选 取 了 CMIP6 气候 模拟 实验 下 的 
ACCESS-CM2 模 式 , 以 及 未 来 试验 场景 SSP2-4.5。 
在 对 2021 一 2030 年 的 ET 模拟 中 展 均 达到 了 0.98 以 
上 。 图 6 图 7 表明 H-Scom 模 型 对 ET 的 季节 循环 以 
及 长 期 变化 趋势 均 有 较 好 的 模拟 能 力 。 因 此 ,H- 
Seom 模 型 在 未 来 也 同样 具有 一 定 的 适用 性 。 


4 讨论 


ET 的 准确 估算 对 于 水 资源 管理 .农业 和 水 文 
气象 研究 以 及 水 资源 收 支 的 确定 至 关 重 要 。 本 


分 区 ”统计 指标 ”修正 前 后 ”最 小 值 ”最 大 值 FHE 
E 区 MAE 原始 2.94 34.00 17.24 
修正 后 1.45 8.48 4.94 
RMSE 原始 3.37 34.58 17.73 
修正 后 1.50 9.59 5.50 
R 原始 0.51 0.90 0.78 
修正 后 0.80 0.96 0.88 
FX MAE 原始 3.43 18.43 10.45 
修正 后 2.62 8.94 4.55 
RMSE 原始 4.01 19.74 11.50 
修正 后 3.40 10.46 5.47 
R 原始 0.60 0.87 0.79 
修正 后 0.81 0.93 0.88 
GK MAE 原始 2.51 15.60 7.91 
修正 后 1.11 3.53 1.94 
RMSE 原始 3.02 15.87 8.51 
修正 后 1.43 3.93 2.37 
R 原始 0.61 0.92 0.82 
修正 后 0.86 0.98 0.90 
HE MAE 原始 9.59 26.94 15.54 
修正 后 1.38 4.00 2.66 
RMSE 原始 10.48 27.17 15.90 
修正 后 1.84 0.79 3.37 
R 原始 0.61 0.79 0.68 
修正 后 0.81 0.96 0.87 


; 尼 为 决定 系数 。 


研究 通过 2 个 研究 区 域 的 128 个 气象 站 点 1961 一 
2020 年 的 逐日 数据 对 H-S 模 型 进行 修正 ,并 进行 未 
来 的 适应 性 分 析 。 

然而 ,尽管 H-Scom 模 型 在 该 研究 区 域 的 精度 有 
所 提高 ,但 对 个 别 精度 较 差 的 区 域 仍 需要 进一步 改 
UE. H H-Scom 模 型 存在 明显 的 空间 变异 性 ,不同 研 
究 区 域 的 时 间 和 空间 尺度 存在 差异 。 本 人 研究 的 时 
空 尺 度 较 大 ,地 理 环境 复杂 ,所 以 在 未 来 的 研究 中 
不 仅 要 考虑 不 同时 间 尺 度 下 的 校正 精度 ,还 需要 对 
不 同 小 区 域 尺 度 的 校正 系数 进行 进一步 的 研究 。 
本 研究 利用 校正 系数 消除 了 风速 .降雨 和 湿度 等 气 
候 因子 对 于 区 域 ET 的 影响 ,但 在 实际 应 用 中 ,这 种 
方法 无 法 直观 地 反映 气象 因子 与 H-Scom 模 型 的 直 
接 关系 ,气候 因子 对 ET 估算 的 影响 还 需 进 一 步 的 
讨论 ,因此 H-Scoma 模 型 不 是 为 了 取代 P-M 模型 , 而 
是 在 只 有 气温 数据 的 情况 下 可 以 更 加 有 效 地 估算 
EToo 
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(c) 2023 年 


(d) 2024 年 
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(e) 2025 年 
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图 6 ACCESS-CM2 模 式 下 2021 一 2030 年 ET, 逐 月 值 与 H-Scomx 模 型 计算 结果 拟 合 
Fig.6 Fitting of monthly ETo values from 2021 to 2030 under ACCESS-CM2 mode and the calculation results of the H-Scorr model 


一 一 H-Scom 模 型 一 一 P-M 模 型 
200 - (a) 2021 年 200 - (b) 2022 年 200 - (c) 2023 年 200 - (d) 2024 年 200 p (e) 2025 年 
€ 150 150 150 € 150 150 
号 100 T a E 100 SI \ El 100 = 100 PN E 100 有 
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246 81012 24681012 2 46 81012 2 4 6 81012 2 4 6 81012 
Att Att 月 份 月 份 月 份 
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图 7 ACCESS-CM2 模 式 下 2021 一 2030 年 ETo 逐 月 值 与 H-Scom 模 型 计算 结果 对 比 
Fig.7 Comparison of monthly ET, values from 2021 to 2030 under ACCESS-CM2 model and the 


calculation results of the H-Scorr model 


来 气候 条 件 下 依然 具有 很 强 的 适用 性 。 因 此 , H- 
5 结论 sum 模型 可 用 于 研究 区 域 的 蒸 散发 预报 ,有 助 于 进 
-RRR DGR. 
(3) 在 长 江 中 下 游 平原 和 西北 黄河 流域 的 具体 
验证 应 用 中 ,其 绝对 误差 结果 表现 参差 不 齐 , 但 整 
体 的 绝对 误差 都 显著 缩小 ,在 西北 黄河 流域 的 验证 
中 效果 优 于 长 江 中 下 游 平原 ,因此 H-Scom 模 型 在 丁 
北 黄河 流域 具有 更 佳 的 效果 。 


(1) 经 修正 后 ,H-Scornr 模 型 具有 较 高 的 精度 ,其 
估算 的 ET, 值 更 接近 于 P-M 模 型 。 长 江 中 下 游 平 原 
MAE FU RMSE 的 平均 值 分 别 下 降 了 6.21 mm 月 "和 
6.38 mm 月 "。 西 北 黄河 流域 MAE FI RMSE 的 平均 
值 分 别 下 降 了 9.26 mm: A! Fil 9.23 mm 月 ,2 个 研 
究 区 域 的 大 部 分 月 份 的 尼 均 大 于 0.9 ,少数 月 份 的 尼 
大 于 0.8, 校 验 后 的 尼 比 校正 前 最 少 提高 1%。H-S 2h (References) 
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Applicability of reference crop evapotranspiration calculation based on 
Hargreaves-Samani regression correction 
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Abstract: The Hargreaves- Samani (H-S) model for calculating reference crop evapotranspiration (ETo) was re- 
gressed and corrected using daily meteorological data from 128 meteorological stations in the Northwest Yellow 
River Basin and the middle and lower reaches of the Yangtze River from 1961, 1981, to 2010. The calculation ac- 
curacy of the H-S upgraded model H-Scorr model was evaluated using site data from 2011 to 2020 and the Pen- 
man-Monteith model was used as the reference. The ACCESS-CM2 model and the future test scenario SSP2-4.5 
are also used to assess the H-Scorr model’ s future adaptability as part of the Sixth International Coupling Model 
Comparison Program (CMIP6) climate simulation experiment. The results show that the mean absolute error 
(MAE) of the four subdomains in the middle and lower reaches of the Yangtze River decreases from 2.58—24.28 
mm to 1.53-7.99 mm per month after correction, and the root mean square error (RMSE) per month decreases 
from 3.22-24.56 mm to 1.96-9.27 mm during the validation period. Monthly, both MAE and RMSE decrease by 
6.21 mm and 6.38 mm, respectively. The MAE of the four subregions of the Northwest Yellow River Basin de- 
creases from 2.51-34.1 mm per month to 1.11-8.94 mm per month, and the RMSE decreases from 3.02-34.58 
mm to 1.43-10.46 mm per month. MAE and RMSE monthly averages decrease by 9.26 mm and 9.23 mm, respec- 
tively. Most months in the two study areas have a coefficient of determination (R°) greater than 0.9, whereas a 
few months have an R greater than 0.8, and the corrected R° value is at least 1% higher than before the correc- 
tion. The R’ values in the CMIP6 climate model’ s future climate scenarios are greater than 0.98, showing strong 
adaptability. Therefore, the H-Scorr model’ s performance in the Northwest Yellow River Basin, as well as the 
middle and lower reaches of the Yangtze River, is improved; it can now better simulate the seasonal cycle and 
long-term trend of ET», as well as improve the accuracy of reference crop evapotranspiration calculations. 

Key words: reference crop evapotranspiration; Hargreaves-Samani model; Penman-Monteith model; model val- 


idation; applicability 


